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ABSTRACT 

Using infrared imaging from the Herschel Space Observatory, observed as part of the Very Nearby 
Galaxies Survey, we investigate the spatially resolved dust properties of the interacting Whirlpool 
galaxy system (NGC 5194 and NGC 5195), on physical scales of ~ 1 kpc. Spectral energy distribution 
modelling of the new infrared images in combination with archival optical, near- through mid-infrared 
images confirms that both galaxies underwent a burst of star formation ^370-480 Myr ago and pro- 
vides spatially resolved maps of the stellar and dust mass surface densities. The resulting average 
dust-to-stellar mass ratios are comparable to other spiral and spheroidal galaxies studied with Her- 
schel, with NGC 5194 at log(M dust /M*) = -2.5 ± 0.2 and NGC 5195 at log(M dust /M*) = -3.5 ± 0.3. 
The dust-to-stellar mass ratio is constant across NGC 5194 suggesting the stellar and dust components 
are coupled. In contrast, the mass ratio increases with radius in NGC 5195 with decreasing stellar 
mass density. Archival mass surface density maps of the neutral and molecular hydrogen gas are also 
folded into our analysis. The gas-to-dust mass ratio, 94 ± 17, is relatively constant across NGC 5194, 
although there is some suggestion that it decreases radially but not significantly above our uncertain- 
ties. Somewhat surprisingly, we find the dust in NGC 5195 is heated by a strong interstellar radiation 
field, over 20 times that of the ISRF in the Milky Way, resulting in relatively high characteristic 
dust temperatures (~ 30 K). This post-starburst galaxy contains a substantial amount of low-density 
molecular gas and displays a gas-to-dust ratio (73 ± 35) similar to spiral galaxies. It is unclear why 
the dust in NGC 5195 is heated to such high temperatures as there is no star formation in the galaxy 
and its active galactic nucleus is 5-10 times less luminous than the one in NGC 5194, which exhibits 
only a modest enhancement in the amplitude of its ISRF. 

Subject headings: galaxies: individual (NGC 5194/M51a, NGC 5195/M51b), galaxies: general, galax- 
ies: fundamental parameters, galaxies: stellar content, infrared: galaxies, tech- 
niques: photometric 
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The fundamental building blocks of a galaxy - its stars, 
gas and dust - are expressed across its spectral energy 
distribution (SED). Hot gas (> 10 6 K) emits brightly 
at very short wavelengths in the X-Ray and ultraviolet 
(UV), while neutral and molecular gas are seen through 
emission lines primarily studied in the submm and radio. 
Young stellar populations dominate the UV and blue vis- 
ible light of a galaxy and ionize gas leading to numerous 
emission lines in the UV, optical and infrared. Evolved 
stars peak at longer wavelengths in the red visible and 
near- infrared regime. At longer wavelengths and lower 
energies, mid-infrared light traces very small (<0.01 /im) 
dust grains and small organic molecules known as poly- 
cyclic aromatic hydrocarbons (PAHs) and far-infrared 
light traces re-processed stellar light emitted by dust 
grains. Multi- wavelength studies provide the synergy 
to compare these components and move toward build- 
ing models that incorporate all the building blocks of a 
galaxy. 

Standard stellar population techniques can tell us how 
much stellar mass is in a galaxy, its metallicity, the 
amount of dust attenuation and can even shed light on 
the history of star formation within a galaxy. Funda- 
mentally, stellar emission models are built from the su- 
perposition of multiple simple stellar populations where 




Fig. 1. — Image of the Whirlpool galaxy system, NGC 5194/5 
(M51a/b) from the Spitzer Space Telescope's IRAC camera at 
3.6 /im (in intensity units of MJysr - 1 ). Measured apertures are 
shown for each galaxy and were defined by a constant stellar mass 
surface density of 50 Mq pc -2 . 

assumptions regarding the IMF and stellar evol ution lead 



to sy stematic but known uncertainties (Conroy et al 



2009). On the other hand, modelling the dust emission 



m a galaxy is not ye t as evolved. Most dust emission 
mode ls or templates (Dal e fc Helou||2002 Galliano et al. 
20111) measure the total dust mass and some dDesert et 
aflfS 



390 



^ Draine fc Li||2007| |Compiegne et al.pOTip the 
relative abundance ol various dust emission components 
(eg. PAHs, very small grains and big grains). Funda- 
mental to these models is the mechanism for dust heat- 
ing, primarily provided by the interstellar radiation field 
(ISRF) produced by stars, although active galactic nuclei 
(AGN) can also play a role in nuclear regions. 

Nearby galaxies offer unique laboratories to explore 
how dust emission is tied to the underlying stellar pop- 
ulations in the context of standard SED modelling tech- 
niques. The post-interacting Whirlpool galaxy system 
was observed as part of the Very Nearby Galaxies Survey 
(VNGS; PI Christine Wilson), a guaranteed-time pho- 
tometric and spectroscopic survey of a diverse range of 
13 nearby galaxies with HerscheVs Photodetector Array 
Camera and Spectrometer (PACS; [Poglitsch et al.||2010[ ) 
and Spectral and P hotometric Imaging Receiver (SPIKE; 
|Griffin et al.|2010| ) instruments. 

This system consists of the diverse galaxy pair of 
NGC5194/M51a, a spiral SABbc type galaxy, and 
NGC 5195/M51b, a dusty spheroid galaxy whose pecu- 
liar dust lanes have gi ven it an SBOp classification (|de| 



Vaucouleurs et al.| l995| ). An encounter between the two 
galaxies 300-500 Myr ago i s inferred from both kinematic 



Salo & Laurikainen 2000) and hydrodynamic modelling 



Dobbs et al. 2010) ol the M51 system. In agreement, 



archival images, suggests a burst of st ar formation oc- 
curre d 390-450 Myr ago in NGC 5195 flTikhonov et al. 
2009). Star formation also peaked around 380-450 Myr in 
the spiral NGC 5194 and has continued until the present 
day. While star formation is ubiquitous and ongoing 
across NGC 5194 /M51a as shown by its - 9 600 HII re- 



gions ([Lee et al. 2011) and optical colors (Tikhonov et 
aL||2009[ ), the spheroid neighbor NGC 5195 shows no ev- 
IcTence of re cent star formation th rough a deficiency of 
ionized gas flThronson et al.| |T991). This is further sup- 
ported by evidence lor dynami cally stable molecula r gas, 
which is unable to form stars. Kohno et al. ( 2002| ) mea- 
sured both CO (1-0) and HCN(l-O) in the galaxy and 
found an HCN-to-CO intensity rati o 5-15 times smaller 
than that seen in starburst regions ( [Heifer fc Blitzp 993; 
Kohno et al.|2002|) an d 5 times smaller than tEa t usually 



found in spiral disks (|Helfer fc Blitzp995| |Kohno et al 
1999), indicative of a lack ol dense molecular cores. 

This paper takes advantage of a large multi- wavelength 
dataset encompassing observations of the stars, dust and 
gas at a common spatial resolution to try to decode the 
star formation history of the galaxy and how it has im- 
pacted the dust and gas properties. Our ultimate goal 
is to present a story of these components in the context 
of the interaction and star formation history of the sys- 
tem. Our observations and the steps involved to match 
the spatial and surface brightness distributions are given 
in Q We apply standard SED fitting techniques to the 
dust and stars separately but do so on a pixel level in or- 
der to map any gradients and spatial dependencies of the 
parameters. The models, our fitting procedure and out- 
put parameter maps are all presented in ^3] We provide 
dust-to-stellar and gas-to-dust mass maps and ID radial 
plots in Qand discuss how the observed dust properties 
relate to the stellar and gas properties of the system. 

2. OBSERVATIONS AND IMAGE PROCESSING 

Our multi- wavelength analysis includes twenty images 
of the Whirlpool system, and thus exploits the great ef- 
forts of a number of large legacy surveys. In this section 
we provide a summary of the archival observations and 
describe our new Herschel observations from the Very 
Nearby Galaxies Survey (VNGS). A number of steps were 
taken to ensure each pixel represents a common mea- 
sure of surface brightness. We have opted to exclude the 
500 /iin image as a compr omise for higher spatial resolu- 
tion. Boselli et al. (2012[) demonstrate that for late- type 



recently constructed color-magnitude diagram of its in- 
dividual stars, measured using Hubble Space Telescope 



spiral galaxies the SF1RE FIR colours are tightly corre- 
lated and the 500 fim band is not necessary to constrain 
a dust model. In addition, because our analysis hinges 
on comparing a large number of observations (4000 pix- 
els in 20 images) to large libraries of stellar and dust 
models, we require satisfactory signal-to- noise levels in 
each pixel to determine parameter likelihoods. In this 
section we describe image processing steps and discuss 
our efforts to accurately quantify uncertainties in post- 
processed images. 

2.1. Adopted Galaxy Parameters 

For ou r analysis we use a co mmon distance of 9.9 ± 
0.7 Mpc (Tikhonov et al. 2009), although a wide range 
of distances (anywhere from 6-10 Mpc) can be found in 
the literature (e.g. 8.4 Mpc from Feldmeier et al. 1997 
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TABLE 1 
Structural Properties 



Galaxy 


RA (J2000) 


DEC (J2000) 


2a 


2b 


P.A. 


Aperture Radius 




(deg) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


(arcsec) 


NGC5194 


202.47065 


+47.19517 


383 


286 


50 


295 


NGC5195 


202.49726 


+47.26531 


285 


225 


90 


115 



is often used). This latest distance is from Hertsprung- 
Russell diagram matching using individual red giant stars 
seen in HST archival images. We also measured our own 
aperture and galactic center for each galaxy based on 
constant stellar mass density rather than surface bright- 
ness. Our apertures, shown in Figure [I] were defined at 
a constant stellar mass of 5OM pc -2 and ensured that 
neither galaxy's aperture overlapped. When total mass 
estimates are measured for each galaxy, it is the sum 
of pixels within this aperture. No aperture corrections 
are applied extending to larger galactic radii because be- 
yond our chosen apertures, the galaxies begin to blend. 
For this reason, mass surface density ratios are more rele- 
vant to our analysis. We also use the apertures to define 
an effective radius, r e ff, in order to measure ID radial 
profiles of mass densities and other derived SED proper- 
ties. Our effective radius is defined as 



7"eff 
fapp 



(i) 



where a is the semi-major axis, b is the semi-minor 
axis, and r^ pp = a 2 + b 2 . The position angle, P.A., is 
measured east from north. These parameters are given 
for reference in Table [I] and shown in Figure [I] imposed 
on a Spitzer /IRAC 3.6 jam image of the galaxies. 

2.2. Far- infrared data 

The VNGS observed the Whirlpool system with Her- 
schel /P ACS at 70 /im and 160 /im and Herschel /SPIRE 
at 250 /im and 350 /mi. Images at each wavelength at 
their native resolution are shown on the left column of 
Figure [2] More details on th e observational setup can be 



found in lBendo et al.| fl2012| ). 

The PA(JS images were processed using both HIPE v5 
and SCANAMORPHOS v^f] To account for changes 
in calibration from HIPE v5 to HIPE v6, the 70 and 
160 /im images were divided by factors of 1.119 and 
1.174. The PACS images were converted to Jy sr -1 from 
Jy pix -1 with platescales of 2" and 4" for the 70 jam and 
160 /im images respectively. 

The SPIRE images were processed using; HIP E and 
BriGAde as described in [Smith et all ( |2012[ ). The 
final SPIRE maps were created using tne naive map- 
per provided in the standard pipeline with pixel sizes 
of 6"and 8"at 250 /im and 350 /im. The FWHM of 
the SPIRE beams for this pix el scale are 18.2" and 
24.5" at 250 /im and 350 /im flSwinyard et al.||2010[ ). 
In addition, the 350 /im data are multiplied by 1.0067 
to update our flux densities to the latest v7 calibra- 
tion product. The SPIRE images required conversion 
from Jy beam -1 to Jy sr -1 where beam is the PSF 
area size of the instrument at the observed wavelengths 

13 http://www2.iap.fr/ users/roussel/herschel/ 



taken for this analysis as 423 and 751 arcsec 2 beam -1 
for the 250 jam and 350 jam images respectively. SPIRE 
images are initially in point source relative spectral re- 
sponse function (RSRF) weighted intensities. To convert 
to extended source RSRF-weighted intensities we divide 
each image by the K^ p terms given in the SPIRE hand- 
book of (1.0113, 1.0087) for the 250 jam and 350 jam im- 
ages. The images remain in RSRF-weighted intensities 
(ie. not converted to monochromatic intensities or color- 
corrected). Model band fluxes are computed by integrat- 
ing the model spectra with each filter's spectral response 
function and thus account for the filter's effective wave- 
length and throughput, as well as the galaxy's varying 
spectral slope. 

The sky background in each image is measured using 10 
pixel wide boxes surrounding the system. For all of our 
Herschel /PACS and i7er.sc/ie//SPIRE images, we sub- 
tracted the median value of the sky region pixels as we 
found values that were on order 1-2 a higher (in the case 
of PACS images) or lower (in the case of our SPIRE im- 
ages) than the standard deviations in the background 
pixels. All images are matched to a com mon resolution 
using kernels created from the PSFs from Ani ano et al 
(2011). They measured azimuthally averaged empirical 
PSFs for the Herschel instruments. The images were 
all convolved to a matched gaussian PSF (FWHM=28") 
slightly higher than the SPIRE 350 jam image to en- 
sure the power in the non-gaussian beam of SPIRE is 
conse rved to the assumed gaussian beam (Aniano et al. 
2011). All images are remapped to a common platescale 
of 10" pix -1 using the WCS tools program REMAFp^j The 
resulting processed images are shown in the middle panel 
of Figure [2] 

2.3. Optical and near/mid-infrared data 

Our optical (BV RI Ha) and mid-infrared (IRAC 3.6, 
4.5, 5.6, 8.0 and MIPS 24 jam) observations come from 
the Spitzer Infrared N earby Galaxies Survey (SINGS; 



Kennicutt et al. 



2003 ). Near-infrared JHK imaging is 
fro m the 2MASS Large Galaxy Atlas survey ( |Jarrett et] 
al. 2003). All images are converted to units of Jy sr -1 



and we account for loss of light due to foreground dust 
attenuation from the Milky Way using the extinction val- 
ues provided by NEE0 using maps from |Schlegel et al.| 
(1998) for all bands up to and including the IRAC 3.6 jam 



image. 

We replace emission from foreground stars (in 10" di- 
ameter apertures) with the median value of emission of 
the pixels contained in its surrounding annulus (4" wide). 
The stars are identified from the 2MASS Point Source 
catalogue. In order to ensure we do not remove any HII 
regions belonging to the galaxy we only remove point 



14 http: / / http: / /tdc- www. harvard.edu/wcstools/ 

15 http://ned.ipac.caltech.edu/ 
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Fig. 2. — Far-infrared images of NGC 5194/5 from Herschel /PACS and SPIRE as part of the VNGS survey from top to bottom at 70 /im, 
160 jUm, 250 /im and 350 fim. The left panel is the original image, while the middle panel shows the image after matching it to a common 
resolution of 28 /; and platescale of 10" pix - 1 . Both images are shown in MJy sr — 1 and the beam area for each image is shown in the bottom- 
right. The right shows the resid ual image bet ween the models and observations, normalized by the observations: (/A,obs — ^A,model)/^A,obs- 
The best-fit model images from |Draine fc Li| ( |2007| > match the images at 160 /im and 250 /im to within ~10%, but underestimate emission 
in some regions of the 70 /im (top) and '6bU /2m (bottom) images by over 20%. 
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Fig. 3. Left: HI mass surface density image of NGC 5194/5 from the THINGS survey ( | Walter et al.|2008|> convolved to m atch the PSF 
(FHWM=28 // ) of our other observations. Middle: the CO (1-0) image from the Nobeyama 45 m telescope (jlvoda et al.|201l|) 



,|201l| ) converted to 



an H2 mass surface density image. Right: Total gas mass surface density image is created by summing the two lert images. All images are 
in units of Mq pc — 2 and are shown on the same intensity scale. Molecular gas dominates the gas emission in the central area of the galaxy 
and along the spiral arms. To convert to total gas mass surface densities, a 1.36 scaling is required to account for helium. 

sources with mid-infrared to optical light ratios compa- 
rable to stars, taking the criterium that the 8 /im surface 
brightness cannot be greater than 1.5 times the i?-band 
surface brightness. If the ratio is higher, Is/Ir > 1.5, the 
point source is an HII region (identified by its bright PAH 
emission at 8 /im) within the galaxy and is not removed. 

Also, we correct for any constant in the sky background 
in each image by placing 10 boxes, each 10 pixels wide 
throughout the sky in each image. If the median in this 
background value is larger than the standard deviation 
across the boxes, we subtracted this constant from the 
image. The narrowband Ha image (not continuum sub- 
tracted) was the only image which required this step. 
Finally, we convolve the images with appropriate kernels 
to our chosen common resolution of 28" assuming gaus- 
sian PSFs for the optical and near-infrared images. The 
FWHMs of the PSFs for the optical images are 2.0", 1.9", 
1.6", 1.9", and 1.8" for the B,V,Ha,R, and I bands re- 
spectively. For the NIR images, the FWHM of the PSFs 
are 3.3", 3.2", and 3.5" for the J, H, and K bands. 

The Spitzer images ar e convolved to a 28 " gaus- 

sian using kernels from |Aniano et al. (2011). The 

FWHMs of the Spitzer images are OO^TSr, 2.11 // , 

2.82" and 6.43" for the IRAC 3.6 /im, 4.5 /im, 5.6 /mi, 

8.0 /im and MIPS 24 /im images respectively. All images 

v/ • -1 



are remapped to a common platescale of 10" pix x using 



the WCS tools program remapj b [ For reference, we have 
included images of the galaxy in each band, before and 
after image processing in the appendix of this paper. 

2.4. Neutral and molecular gas 

Imaging of the neutral ga s is provided by the HI Nearby 
Galaxy Survey (THINGS; [Walter et al.|[2QQ8| ). Already 
in units of M pc -2 , we match plate scales and convolve 
the beam assumed to be a gaussian of FWHM=6" to a 
28" beam to match our other observations. Molecular gas 
observ ations at CO (1-0 ) from the Nobeyama 45 m tele- 
scope (Koda et al. 2011[ ) are also included in the analysis. 
The CO (1-0) image, originally in antenna temperature 



http: / / http: / /tdc- www. harvard.edu/wcstools/ 



units of K kins" 1 , is divided by 0.4 to convert to main 
beam temperature units. We then use a constant conver- 
sion factor of X C p = 2.0 x 10 20 mol cm -2 (K km s -1 ) -1 
( Strong et al.||l988l) which correspo nds to 3.2 M pc -2 
(K km s" 1 )"" ( [Wilson et al.||2QQ9| | to convert the CO 
intensities to an H2 gas mass surface density. This im- 
age is then convolved from a 22" beam to a 28" beam 
to match all of other observations as well as changing 
the platescale to 10" pix -1 . Both images of the gas mass 
surface densities are shown in Figure [3] The sum of both 
the neutral and molecular gas is shown on the right. All 
images are multiplied by 1.36 to account for helium and 
thus represent the total, not just hydrogen, gas in the 
galaxy. We use the total gas mass surface density in our 
comparisons to the stellar and dust mass surface densi- 
ties. 

2.5. Uncertainties and qualifying signal-to-noise criteria 

We generally considered two types of uncertainties in 
each pixel value of each image: the error due to the un- 
certainty in the sky background and also the calibration 
uncertainty. The uncertainty in the sky background is 
measured in each image after background subtraction 
and convolving to lower spatial resolution by selecting 
ten 100" wide boxes throughout the image. The stan- 
dard deviation in the median value of the pixels in each 
box leads to our uncertainty in the sky background. Cal- 
ibration uncertainty is taken from the various instrumen- 
tal handbooks or survey description papers. Calibration 
uncertainties in the optical data are assumed to be 5%, 
except for a 10% uncertainty in the narrowband Ha im- 
age. Both 2MASS and Spitzer images have calibration 
uncertainties on the order of 3%. Uncertainties for Her- 
schel observations are 3% and 5% at 70 /im and 160 /im 
and 7% at 250 /im and 350 /im. In addition, error maps 
of the Herschel images were also included in the analy- 
sis. These images were processed to a common spatial 
scale as the images above by taking the square root of 
the convolution of the square of the error image with the 
square of the kernel. Furthermore, the error images were 
matched to the 10"/pix platescale. Uncertainties in each 
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Fig. 4. — We plot four example SED fits for single representative pixels from the nuclear region (top-left), inner arm (top-right) and 
interarm (bottom-left) regions of NGC 5194 and in the bottom-right, a pixel from the center of the post-starburst NGC 5195. Model 
band fluxes of each best-fit SED are shown in open triangles. The purple curve is the best fit stellar emission model (PEGASE.2, [Fic^l 
|fc Rocca-Volm erange 1997) to the B-K images (Ha excluded) and the red curve is the best fit dust emission model from DL07 to the 
5. t>-35U jiim images, we separate the PDR component (dotted orange line) from the dominant ISM component (dashed line) for each DL07 
model. High PDR fractions are found on the spiral arm of NGC 5194 (top-right) and in the central region of NGC 5195 (bottom-right). 
For comparison, we also show the best fitting modified greybody (dashed green curve) for both examples and give their characteristic 
temperatures and dust mass. Note that the dust and stellar SEDs have not been added together and neither SED was constrained by the 
IRAC images at 3.6 /im and 4.5/im. 



pixel were added in quadrature to produce a combined 
error map at each bandpass. 

We want to include as many pixels in our analysis as 
possible and at the same time exclude regions in which 
the majority of pixels have low signal-to-noise (S/N) lev- 
els. A fit to the SED of a pixel was performed if at least 
5 bands had S/N> 3 in the optical/near-infrared data 
from B through K-band. We also required S/N> 3 in 
at least 3 bands for the infrared data from 5.6 jam to 
350 /im. In practice this selection was dominated by the 
relatively less sensitive infrared data; every pixel which 
satisfied the infrared criteria had S/N> 3 in all optical 
bands. 



3. SED LIBRARIES AND FITTING RESULTS 
Although models exist (eg 



Cunha et al. 



Devriendt et al. 



1999 



2008[ |Noll et al.||2(JU9| jJ^p escu et al.||2(jil| [ 



da 



which describe both stellar and dust emission together, 
they tend to be either computationally expensive (re- 
quiring many fitted parameters and slow codes comput- 
ing the uncertainties in these parameters) or lack the 
most up-to-date SEDs. Unlike today's commo n stellar 



Rocca-Volmerange 


1997), GALEXEV Bruzual & Char- 


lot||2UU3|)j and FSJ 


ts (IConroy et al. 2UU9|) - these sell- 



consistent stellar and dust emission codes are not exten- 
sively tested and their use is not widespread. For this 
reason, we opt to fit the stellar SED and dust SED com- 
ponents separatel y. This approach is similar to that of 



Noll et al.| ([2009) but by considering both components 
separately we can fit the SED at a much faster speed; 
this is primarily out of necessity since we are dealing with 
a very large dataset (20 images with 4000 pixels each). 
The downfalls of this approach are first that energy is 
not necessarily conserved from the optical to IR (such 
that the re-processed infrared emission is linked to the 
ionizing radiation provided by th e stellar populations) , 
something that both MAGPHYS (Ida Cunha et al.|2008 ) 
and CIGALE flNoll et al.1|2009D do achieve and second 
the near- to mid- infrared regime from ~ 1 — 5 /im where 
both dust and stellar emission is significant cannot be 
used to constrain the fit. However, because both pre- 
main-sequence and post-main-sequence stellar emission 
is not well calibrated in the SPS models at these wave- 
lengths, nor is PAH emission, it is not ideal to fit this 
regime in any case. 

In our method, we first fit our B through K band pho- 
tometry to stellar emiss ion models using standard SPS 
codes as we describe in £3.1 Using the predicted stellar 
emission from these fits, we are able to subtract stel- 
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Fig. 5. — Stellar mass surface density in Mq pc 2 is shown in the left panel and the normalized uncertainty in the pixel mass determination 
is shown in the right panel as M rms /M av g- 



lar emission from the mid-infrar ed im ages at 5.6 /im and 
8.0 /im. As we will describe in §3.2[ we then fit the 5.6 
through 350 /im images to the mid- through f ar-infrared 
dust emission models from Draine & Li ( 2007| ) (hereafter 
DL07), as well as fitting a single temperature modified 
blackbody model. 

3.1. Stellar SEDs 

A number of codes exist for a user to generate a model 
galaxy SED. The choice of parameters and assumptions 
that are entered into the codes can lead to systematic 
differences in the output parameters. Stellar masses are 
most robust and can be predicted to w ithin 0.2 dex 



(Muzzin et al. 2009 Conroy et al. 2010) and often to 



higher precision it near- infrared photometric bands are 
included, although masses will vary systematically de- 
pending on the choice of initial mass function (IMF). 
Other paramete rs tend to have larg er dispersions across 
different codes ( |Muzzin et al. 2009). If uncertainties in 
the inputs (sucn as evolutionary isochrones, stellar li- 
braries, binary fraction) of the code are conside red, then 
mass es are only good to within ~0.3dex ( |Conroy et al.| 
2010|). For a more intensive discussion on systematic 



differences and known uncerta inties of SPS models , the 
reader is referred to paper s by Muzzin et al. ( 2009| ) and 
IConroy et aL| ( [20091 [2Q10| ) . 

In this study we take a more empirical approach and 
try to compare the best fit model SED parameters to in- 
formation from other observations and independent stud- 
For example, we have some prior knowledge re- 



les. 



garding the star forming history of the system based on 
isochrones (which have some of the same uncertainties as 
those found i n SPS codes) fit to ind ividuals stars, most 
recently from |Tikhonov et al. (2009). They indicate an 



enhancement in star formation around 400 Myr ago in 
both galaxies. Independently, previous dynamical stud- 
ies of the system provide a consistent temporal view of 
the system. The motions of gas and stars in the galaxies 
suggest at least one encounter between the two gal axies 
300-500 Myr ago from tw o independent studies from |Salo| 



30-500 Myr ago from tw o independent studies from Salo 
LaurikainenT ( |2000| ) and |Dobbs et al.| ( |2010| ). This prior 
lowledge of the star formation history could be used to 



knowledge of the star formation history - 
constrain our fitting, but instead we have opted to keep 
the star formation history (SFH), stellar population ages 
and dust attenuation as free parameters to see how they 
compare to these independent studies. 

Our stellar libr aries are generated by the exten sively 
tested PEGASE.2 ( |Fioc fc Rocca-Volmerange||1997[) code 
and are an extension of work done in [plaze brook e"tl 



alj fl2004| ) an d |Mentuch et al.| ( 2009 ) and include four 
metalTiclties: Z= {0.004, 0.008, 0112^.05}. We assume 
no evolution in the metallicity and no gas infall, although 
this is likely not the best explanation for NGC5195 
in particular since it appears to be accreting mate- 
rial from NGC 5194. Our adopted IMF is that of 
Kroupa (2001). The SFH is modelled as a sum of 



a continuous star forming population described by a 
tau exponential and an additional burst of star forma- 
tion whose amplitude and age are both free parame- 
ters. The SFH of the continuous population is mod- 
elled as an exponentially decreasing function of time 
(exp -£ / r ), parameterized by an e-folding timescale of 
r = {100, 500, 1000, 5000, 500000} Myr where 100 Myr 
is a short burst of star formation or roughly a sim- 
ple stellar population and the largest timescale of 500 
Gyr represents a constant star formation history. The 
burst component is modelled as a tau-exponential with 
r = 100 Myr. The time since the star formation 
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R eff ( k P°) R eff ( k P°) Reff ( k P C ) 

Fig. 6. — Radial profiles of SED derived star formation history parameters for NGC 5194 in purple dots and NGC 5195 in cyan triangles. 
Left: The age of the old stellar population is logf* = 9.9 db 0.2. Middle: Age of the burst component. A secondary star formation 
burst period happened 300-600 Myr ago. Right: The mass fraction of the secondary burst component relative to the older primary stellar 
population is ~ 5 — 15% throughout the galaxy. Note that the deviation from pixel to pixel tends to be less than the uncertainty in the 
parameters indicating that at the spatial scale of our analysis (28 /; ), the stars are well mixed. 



burst is allowed to be an additional parameter and 
was given the largest range and sampling at steps of 
50 Myr from 50 to 1000 Myr and then steps of 250 Myr 
from 1 Gyr to 5 Gyr and then steps of 500 Myr from 
5 Gyr to 13 Gyr. The age of the galaxy could be 
anywhere between 1 Gyr and 13 Gyr with fairly high 
temporal sampling as well. The mass of the burst 
compared to the mass of the continuous star form- 
ing population (CSP) ranged as log (Mburst A^csp) = 
{-3, -2, -1, -0.75, -0.5, -0.25, -0.1, 0, 0.1}. Thus, the 
burst SSP could be anywhere from 0.001 to 1.25 times 
the mass of the primary CSP. 

As we did not include UV observations in this analy- 
sis, the nuances of the 2175 A bump did not have to be 
considered in the choice of a dust attenuation law. We 
opted to u se the Milky Way dust attenuation curve from 
Pei (1992) and only considered a screen typ e geometry 
of t he dus t. T his law, compared with those of |Calzetti et| 



al 



_ fll994| ) and |Charlot fc Fall] fl2000p for example, more 
accurately describes the dusT attenuation of stars in the 
Milky Way ( |Indebetouw et al.|2005[) ob served in the near 
and mid-infrared ( jlVlentuc h et al. 2010). C urves t hat rep 



resent star b urst galaxies (Calzetti et al. [1994 Chariot 



fc Fall||2000|) incorrectly model the dust attenuation be 
cause they ha 



that is 



ave not accounted for non-stel lar emission 
uite large in star forming regions ( |Mentuch et| 
al. 2010| ) whi ch becomes significant at wavelengths be- 
yond ~1 /im ( |Mentuch et al. 2009 ) . The amplitude of 
the dust attenuation was an additional free parameter 
ranging from Ay =0 — 2 mag in steps of 0.25 mag. 

We place the modelled galaxies at a distance of 9.9 Mpc 
to match the system's distance and integrate the syn- 
thetic SEDs across each observed band's filter response 
curve to get theoretical band fluxes. A least-squares com- 
parison between observed and modelled band fluxes in 
each pixel selects the best fit model from our suite of stel- 
lar models to our BVIRJHK maps. The stellar mass 
in each pixel is inferred from the least-squares normal- 
ization of the model to observed band fluxes. 

Following standard Monte Carlo error simulations, our 
errors are derived for individual parameters by running 
our fitting routine one hundred times, each time allow- 



ing the pixel values in each band to sample the gaussian 
probability distribution of each measurement. The un- 
certainty, or 1 a spread in the gaussian distribution, ac- 
counts for calibration uncertainty, sky backgro und limi- 
tations and other uncertainties as described in £12.51 The 



mean and standard deviation of the best fit parameters 
from all trials lead to our best estimate parameter value 
and its uncertainty in each pixel. 

Figure[4] shows four examples of fits to individual pix- 
els in different regions of NGC 5194/5. The purple curve 
shows the best-fit SED to the B-K pixels. The large de- 
crease in amplitude from the nuclear region of NGC 5194 
(top-left) to the arm regions (top-right and bottom-left) 
is due to a decrease in stellar mass density. Overall we 
find that the models fit the observations to better than 
10% and are consistent with the observations within our 
estimated uncertainties. In an appendix, we show qual- 
ity of fit images for the optical, near- and mid-infrared 
images. Consistently we have very good fits of the stel- 
lar emission, with X? e duced ~ 1 f° r mos t of NGC 5195 
and Educed ~ 0.5 - 1 for NGC 5194. Only the Ha im- 
age (which we don't use to constrain the SPS model) 
is slightly mismatched. We have simplified the SFH in 
the model to consist of just two components. But we 
know that for NGC 5194, it consists of at least two con- 
tinually star forming episodes: one old stellar population 
and then one that has been forming stars since at least 
its encounter with NGC 5195 and likely before given its 
gas fraction and bulge-disk ratio. But our simplified SFH 
only consists of one continually star forming population 
and then a burst population. Because the current star 
formation is not being accounted for properly, neither is 
its associated emission in the Ha narrowband image. 

Maps of the stellar mass surface density and fractional 
uncertainties (M*, rms/M*, avg) are shown in Figure [5] 
On average, stellar masses in each pixel can be found to 
a precision of 0.1 dex based on our Monte Carlo simula- 
tions. The distribution of stellar mass radially decreases 
from the nuclear regions of both galaxies and there is a 
slight enhancement of stellar mass along the spiral arms 
in NGC 5194. 
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Fig. 7. — The left panel displays a map of metallicity across the NGC 5194/5195 system showing increased metallicity in the centers of 
both galaxies, relative to larger radial distances. The stellar mass surface density of the system is shown in contours. Uncertainties on this 
parameter tend to be quite high as can be better seen in the ID radial plot of metallicity for NGC 5194 shown on the right. Grey points 
are for individual 30 /; pixels and the purple points are the pixels binned radially The metallicity tends to decrease with radial distance 
consistent with results from Moustakas et aL] ( |2010| . Their abundance measurements are calibrated using two methods. He re we show 
best-fit lines to the radial abundance profile cTerived usin g two different calibration methods: the Kobulnicky & Kewley (20(M| calibration 
method in dashed red and the |Pilyugin fc Th uan (2005) calibration method in dotted blue. 



Star formation history parameters: In Figure [6] we give 
radial plots of parameters describing the star formation 
history (SFH) of the system. The profiles go out to 
a radial distance of 295" for NGC 5194 and 115" for 
NGC 5195. Beyond these radii the two systems signif- 
icantly overlap. The left panel shows the age of the old 
stellar population. There is no significant radial depen- 
dence and both systems have similar old stellar popu- 
lation ages of log t*[yr] = 9.9 ± 0.2. Thus, the age is 
not well constrained for the system, ranging from ~7- 
lOGyr. Somewhat reassuringly, the SED fits reveal, as 
we show in the middle panel of Figure [6j that there was 
an increased episode of star formation across the sys- 
tem. For the inner radius of 295" of NGC 5194 we find 
a burst occurred 3.8 ± 0.9 x 10 8 yrs ago and for the in- 
ner 115" of NGC 5195 the average time since the burst 
is 4.2 ± 0.8 x 10 8 yrs ago, consistent wit h th e dynami- 
cal mode ls of |Salo fc Laurikainen| ( 2000 ) and Dobbs et 
|al.| pOlO), and the color-magnitude diag rams of individ- 



ual red giant branch stars seen by HST (Tikhonov et al. 
2009). Our analysis suggests the fraction ol mass, rela- 
tive to the primary population, formed in the burst in 
NGC 5194 is 0.12 ± 0.05 and in NGC 5195 is 0.10 ± 0.05. 
A radial profile of the burst fraction for both galaxies is 
shown in the right panel of Figure [6] and again, no spatial 
dependence is found. 

Metallicity: While the SFH parameters are generally well 
mixed in the galaxy at the spatial resolution of this work 
(~ lkpc), the metallicity of the stars shows a trend of 
decreasing metallicity values with radius. In Figure [7| 
we show a map of the metallicity, where the metallicity 
transitions from being super solar in the nuclear regions 
of both galaxies (topping out at Z=0.05, our maximum 
fitted metallicity value), to being near solar (Z=0.02) 
on the outer galaxy edges. The uncertainty from our 
Monte Carlo simulations is quite large because of the 
poor sampling of this parameter available by PEGASE.2. 



The right plot shows the ID radial profile of the metal- 
licity for NGC 5194. Our results are consistent, but not 



as precise, w ith r adial metallicity profiles from |Bresolin 
et al.J Q2004| ) and |Moustakas et al. Q2010 ) who also find 
metallicity decreases with radius in NGC 5194. Most 
late-t ype galaxies exhibit de creasing metallicity gradi- 
ents (M oustakas et aL]|2010[ ), whose slopes depend on 
the calibration method to convert line strengths to abun- 
dances and d iffer f rom galaxy to gala x y. Bo th |Bre solin 



et al. (2004) and Moustakas et al. (2010) measured 
the metallicity radial profile in NGC 5194 using abun- 
dances derived from optical spectroscopy of HII re gions. 
We convert t he be st-fit abundance profiles from |Mous- 
takas et al.| ( |2010[ ) to a metallicity radial profile (as- 
suming Zq = U.Uz) using the value for sol ar metallic- 
ity (12 + \ ofr{Q/H)n = 8.69 ± 0.0 5) from |Asplund et 



al. Q2009 ). Moustakas et al. ( 2010[ ) use two calibration 
methods to convert line strength to chemical abundances. 
We plot both calibration methods in Figure [7[ The red 
dashed line has b een converted following the |Kobulnicky| 
|fc Kewley] ( |2004[ ) c alibration method and the blue dotted 
line following the Pilyugin & Thuan (2005) calibration 
method. Our values lie somewhere n the middle, but 
the same trend in decreasing metallicity is seen. System- 
atic offsets result from our assumed solar metallicity as 
we ll as the calibration m ethod. The reader is referred 
to |Moustakas et al. (2010) for more precise metallicity 
results! 

Dust attenuation: SPS models also include a parame- 
ter that quantifies the intrinsic dust attenuation. Fig- 
ure [8] shows a radial plot of the dust attenuation for each 
galaxy. There is enhanced dust attenuation in the cen- 
tral regions of both galaxies which decreases with radius 
from 0.5 mag of visual dust attenuation in the center 
down to about 0.2 mag. For comparison we plot the ra- 
dial profile of d ust extinction for NGC 5194 measured in 
the near-UV by|Muhoz-Mateos et al.|([2009). We convert 
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Fig. 8. — Radial profiles of visible dust attenuation, A v in units 
of mags shown for NGC 5194 (purple dots) and NGC 5195 (cyan 
trianges). A v is a fitted parameter in our SED stellar popula- 
tion fitting analysis. For compariso n, the orange squares s how the 
measured visual extinction from Muhoz-Mateos et al. 2009 inferred 
from their measured near-UV attenuation. The dashed orange line 
show the systematic uncertainty due to the 2175 A bump. 

from A(NUV) to A(V) us ing a facto r of 3 based on the 
dust extinction models of |Pei| ( |1992| ) and |Draine| (|2QQ3 ) 
but caution that this attenuation ratio is sensitive to the 
UV bump at 2175 A and can range anywhere from ~ 2.4 
to 3.3. The dashed orange lines on Figure [8] demonstrate 
the systematic uncertainty due to this conversion factor. 
Error bars show the precisio n in the measurement from 



|Muhoz-Mateos et al. (|2QQ9[) which is more precise than 
ours and generally higher. Dust attenuation is modestly 
underestimated likely because of exclusion of UV and U - 
band data in our analysis. 

3.2. Dust SEDs 

We use the modelled stellar maps of the galaxy to sub- 
tract the contribution of stellar emission from our mid- 
infrared images of the IRAC 5.6 /im and 8.0 jam bands. 
These resulting images, which probe emission related 
to polycyclic aromatic hydrocarbons (PAHs), and our 
other infrared (24-350 /im) images are fit using the two- 
component emission model put forth by ||Draine &; Li 
( [2007| ) (hereafter DL07) as well as a conventional single 
temperature modified blackbody component. 

3.2.1. Draine and Li 2007 dust models 

DL07 model the mid- through far-infrared emission due 
to PAHs, very small dust grains and larger dust grains 
expressed as a function of the underlying interstellar ra- 
diation field (ISRF). The dust is assumed to consist of a 
mixture of carbonaceous grains and amorphous silicate 
grains, with size distributions that are consistent with 
the observed wavelen gth-dependent dust attenuati on in 
the local Milky Way (We ingartner fc DrainepOOlj ) , and 
allows for varying PAH abundances. The infrared emis- 
sion is a function of the illuminating radiation field which 
in this model is assumed to come from two components. 
Most of the emission comes from diffuse dust heated by 
the ISRF. The bulk of the dust mass in a galaxy is in 



this component. A second component represents emis- 
sion from photodissociation regions (PDRs) where mas- 
sive star formation creates much stronger radiation fields, 
leading to relatively higher emission at shorter infrared 
wavelengths. The diffuse ISRF is scaled to match the 
spectrum of the ISRF of the MW. The amplitude of the 
diffuse ISRF is a fitted parameter, U m i n , normalized to 
the amplitude of the ISRF of the MW. The radiation 
field of the PDR component consists of a range of ener- 
gies from U m in to U ma x, where U ma x can be an additional 
fitted parameter and U m in is the ISRF of th e diffuse com- 
ponen t. We opt to fix Umax = 10 6 Umw as |Draine et al.| 
(2007) showed that it was not a highly sensitive parame- 



ter and that most galaxies from the SINGS survey were 
fit well by this fixed value parameter. 

Our model library consists of MW type dust mod- 
els from DL07 and varies in U m in from 0.1 — 25 x Umw- 
The relative contribution of the secondary PDR compo- 
nent is quantified by the parameter, 7, and ranges from 
0.001 to 1 times the primary ISM component. An ad- 
ditional parameter quantifies the mass fractional contri- 
bution of PAH emission, <2pah, i n the mid-infrared re- 
gion of the SED with following range of values: qpAH = 
{0.47,1.12,1.77,2.5,3.19,3.9,4.58}%. Finally, we con- 
vert the SEDs from units of per hydrogen nucleon to 
units per dust mass assuming Mdust/Mn = 0.01. 

As with the stellar SEDs, model band fluxes are gener- 
ated by integrating the PACS and SPIRE filter response 
functions over the model SED. The SPIRE filter curves 
were multiplied by a normalized A 2 function to convert 
the response curves from those defined for point sources 
to those of extended source emission. A least-squares 
fit is done between the observed bands from 5.6 /im to 
350 /im to obtain the best fit parameters. Again, average 
best fit parameter values and uncertainties are measured 
from 100 Monte Carlo simulations of the probability dis- 
tribution of the photometry. The dust mass in each pixel 
is inferred from the least-squares normalization of the 
model to observed band fluxes. 

Overall, the fits reproduce the 5.6-350 /im emission 
fairly well, leading to X? e duced = 1 — 2 over much of 
NGC 5195 and the central 200" aperture of NGC 5194. 
Outside of r=100 // , the quality of fits degrades, resulting 
in X 2 e duced > 10- This can be seen in the disagreement 
between models and observations at 70 /im and 350 /im in 
Figure [2] The right panel shows the quality of fit of the 
model, defined by the following equation: 



^A,obs — I\ ,model 



*A,ob 



(2) 



In general we find that the models underestimate the 
emission at 70 and 350 /im by as much as 20% but at all 
other wavelengths the photometry and models match to 
better than 5% over most of the galaxy. 

The dust mass surface density map is shown in Fig- 
ure [9] along with the fractional uncertainty in deriving 
the parameter. The dust mass is fairly well constrained 
by the models with the parameter having less than 10% 
uncertainty over most of the galaxy. The dust mass is 
better constrained by the dust models than the stellar 
mass is by SPS models because of the smaller variation 
in the mass-to-light ratio of dust emission. Unlike the 
stellar mass map (Fig.|5|, there is more spatial structure 
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Fig. 9. — Maps of the dust mass surface density in Mq pc -2 (left) with the stellar mass surface density from Figure [E] shown in contours 
at levels of 50, 100, 200, 500, 1000, 1500, 2000, 2500 and 3000 M©pc -2 . The normalized uncertainty determined from our Monte Carlo 
simulations in the pixel mass determination is shown in the right panel as M rms /M av g- 



in the dust mass which is particularly enhanced along the 
spiral arms, similar to the mass distribution seen in the 
molecular gas in Figure [3] NGC 5195 does not show a 
radially decreasing profile as seen in its stellar mass map 
but still does contain a fair amount of dust in its vicin- 
ity and there is some indication of central concentration. 
If we consider the higher spatial resolution images from 
70 /im-350 /mi in Figure [2| we observe that the thermal 
dust emission is very peaked, which could suggest that 
the dust is concentrated. However, the main reason for 
this is not that the dust mass is higher, but rather the 
illuminating radiation is much more intense and concen- 
trated, leading to the peaked infrared emission seen in 
the infrared images. 

The left panel in Figure [lO] shows the U m i n parame- 
ter map of the amplitude of the diffuse ISRF compo- 
nent. The ISRF in NGC 5195 is extremely high, up to 
25 times that of the ISRF in the Milky Way. Over most of 
NGC 5194 the ISRF is very close to the Milky Way value 
(~2-3 Umw), but increases by a factor of 7-15 Umw along 
the spiral arms. The nuclear region is also higher by a 
factor of 10 Umw- Our formal fitting uncertainties in this 
map and the one in the middle of the PDR parameter 
are quite small (about 5%) indicating the infrared data 
was quite sensitive to both parameters. 

The middle panel shows the fractional contribution of 
radiation associated with PDR regions, that is regions 
with intensities as high as 10 6 times that of the MW. 
The radiation coming from PDR regions in NGC 5194 
is concentrated along the spiral arms but not the nu- 
clear region. About 2% of the emission in NGC 5195 can 
be associated with this type of emission, even though 
there is no massive star formation providing the radia- 
tion. The bright emission seen in the infrared images 



comes from intense radiation fields, not enhanced dust 
mass. The high dust heating found in NGC 5195 pro- 
vides strong support for the dust to be located within 
the galaxy rather than in front of it. But the question 
then remains as to what is heating the dust? The stellar 
density is higher in NGC 5195, but not that much higher 
than in the bulge of NGC 5194, so there is another m ech- 
anis m at play. We leave this for later discussion in £4.2 



and §473[ 

Finally, our SED fits constrain the PAH emission and 
is shown in the right panel of Figure [lO] The PAH frac- 
tion is the maximum value allowed by the models of 
4.58% (and thus could feasibly be higher) across most 
of the system. The fraction remains high along the spi- 
ral arms of NGC 5194 extending out to NGC 5195 but 
then decreases across NGC 5195 to about half the value 
3.35 ± 1.07 ) where a compact radio source 



van 



der Hulst et a l.|1988| ) attributed to an AGN exists. There 
is also a small depression in the abundance in the center 
of NGC 5194, which is a low - luminosity Seyfert 2 nucleus 
dTerashima fc Wilson||2001 ). 



3.2.2. Modified black body models 

A modified black body model was also fit to the 70 
though 350 fim images as: 



n v M d B(y,T d ) 
D 2 



(3) 



where M d is the dust mass with a temperature of T d 
at a distance, D. The emission is expressed as the 
Planck funct ion, B(v,Td), modified by n v 



oc v H 



Boselli 



et al. ( 2012[ ) indicate that for at least late- type, metal- 



rich galaxies a value of f3 = 2 appropriately models the 
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Fig. 10. — Left: Dust parameter map of the amplitude of the ISRF of the diffuse dust component, U m i n . Stellar mass density is shown in 
dashed contours in all three panels at levels of 50, 100, 200, 500, 1000, 1500, 2000, 2500 and 3000 M pc~ 2 . The ISRF is enhanced along 
the spiral arms by a factor of ~ 7 — 15 in NGC 5194. Surprisingly, U m i n is highest in the post-starburst galaxy, NGC 5195, and is up to 25 
times the amplitude of the ISRF of the Milky Way. Middle: Parameter map of the fractional contribution of emission from PDR regions to 
the total ISRF. The PDR fraction is highest in NGC 5195 even though there is no star formation. Right: Mass fractional contribution of 
PAH molecules to the total dust mass. Over most of the system, the abundance is 4.58%, the highest allowed in the models. Both nuclear 
regions show small depressions where both harbor low-luminosity AGN. 
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Fig. 11. — Dust mass surface density (left) and temperature 
(right) from single temperature modified blackbody models. Just 
as with the ISRF parameter, U m i n from DL07, the dust tempera- 
ture from modified black body fits is largest in the post-starburst 
galaxy NGC 5195. Temperatures are slightly higher in the nuclear 
region and along the spiral arms. 



SPIRE infrared colors, but caution that ft = 1.5 appears 
to be more appropriate for less metal-rich, lower surface 
brightness lat e-type galaxies in t he Herschel Reference 
Survey (HRS; |Boselli et al. 2010). As a test, we verify 
that a $ = 2 value is sufficient. We performed simi- 
lar modified blackbody fits allowing f3 to vary from 1 to 
3, but this time fit at lower resolution and include our 
500 fim SPIRE image. We found that /3 = 2.0 ± 0.6 over 
the whole system. 

We normalize for the dust mass at 250 /am with 
^(A=250/im) = 0.398 m 2 kg _1 , using the opacity func- 
tion for the Milky Way from |Draine| (|2003), updated 
from their dust model presented in Weingartner & Draine 
( |2001| ) and consistent with the opacity model from the 
D" 



although i t should be noted that recent results from Her- 
schel (e.g. |Bendo et al. 120 10a[ [2012] |Boquien et al.|2011 ) 
suggest such a simplistic model is not appropriate since 
the dust emitting at < 160 /im is not the same thermal 
component as the dust emitting at > 250 /im in nearby 
spiral galaxies. 

As above, band fluxes were generated by integrating 
the model SED over the filter response curves appro- 
priate for extended source emission. The normalization 
between the model and observed band fluxes gave the 
dust mass. Even though our models were normalized at 
250 /im, the fits themselves were not particularly biased 
towards this band since our least-squared normalization 
yields the dust mass not a single band flux. The infrared 
observations were sufficiently modelled by the single tem- 
perature model resulting in X? e duced ~ 1 over the entire 
system, even though other galaxies are not so easily de- 
scribed. In fact, single temperature models fit all the 
images from 70-350 /im as well as the DL07 models if fit 
to the same band fluxes. 

The resulting dust mass distribution, seen in Figure [TTj 
is fairly similar to the one derived using the DL07 mocL 
els from our Figure [9| The temperature distribution also 
qualitatively matches the ISRF amplitude map from the 
DL07 models from our Figure [10] Although the qual- 
itative structure is similar, quantitative pixel-by-pixel 
comparison of the dust mass shows that systematically 
modified blackbody masses overestimate the dust mass 
relatively to the DL07 models at high ISRF values. In 
the left panel of Figure Il2l the ratio between derived 



)L07 models. We assumed a single temperature model, 



dust masses from the DL07 models and modified black- 
body models are shown for each galaxy as a function of 
the amplitude of the ISRF. At higher ISRF, the light-to- 
mass ratio of the dust is higher. In the central region of 
NGC 5195 where the ISRF peaks, the dust mass is over- 
estimated by about a factor of 3. Without accounting 
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Fig. 12. — Left: DL07 dust masses relative to single dust tem- 
perature dust masses. Right: Dust temperature from modified 
blackbody fits is compared to the DL07 ISRF parameter, U m i n . 
As expected, the two parameters are logarithmically related. 

for this variation in the light-to- mass ratio, the modified 
blackbody masses overestimate the dust mass relative to 
the DL07 dust masses. This is analogous to the optical 
regime in which young stellar populations exhibit much 
larger light-to-mass ratios (particularly in the blue) than 
evolved populations. We also show the relation between 
dust temperatures and the DL07 U m i n from the different 
models in the right panel. As expected, the two are log- 
arithmically related, although a few values diverge from 
the trend in NGC 5195. 

4. DISCUSSION 

In the previous section, we presented the output pa- 
rameter maps resulting from the average best-fit parame- 
ters from 100 Monte Carlo simulations of the stellar emis- 
sion and dust emission fitting routine. Here we compare 
some of the parameters and discuss their implications. 
First we provide mass comparisons and ID radial profiles 
of the stella r, ga s and dust mass surface density for both 



galaxies in £4.1 We then turn our discussion to the re- 
lation of interstellar heating of the dust in the context of 
other properties of the galaxies in £4.2 We focus on the 
peculiar properties of the post-starburst NGC 5195 with 
its high amplitude ISRF, ~30K c entra l dust tempera- 
ture and lack of star formation. In §4.3| we synthesize its 
properties and discuss why it is no longer forming stars 
despite containing molecular gas. 

4.1. Mass comparisons 
4.1.1. Radial projections and total masses 

In Figure [l3j we plot the radial mass surface density 
profiles of the stars, gas and dust in black points, blue tri- 
angles and red squares respectively for NGC 5194 in the 
left panel and NGC 5195 in the right panel. The effective 
radius, R e ^ , is determined using the structural parame- 
ters listed in Table[l] The apertures of both galaxies 
were chosen at a constant stellar mass surface density of 
5OM pc -2 . This threshold ensured that the apertures 
were not overlapping. 

In the spiral galaxy, NGC 5194, the profiles of the stars, 
gas and dust all show similar slopes, but in NGC 5195, 
the slope of the stars is slightly steeper and there is no 
significant slope in the gas and dust distribution. As can 
be seen in the gas and dust mass surface density plots in 
Figures [3] and [9j the gas and dust in NGC 5195 appear to 
be connected to the spiral arm of NGC 5194, rather than 
centrally concentrated and coupled to the stellar mass 





M* (10 10 M ) 


M gas (10 10 M ) 


M dust (10 8 M ) 


NGC 5194 


4.7 ±0.1 


1.10 ±0.03 


1.19 ±0.01 


NGC 5195 


2.5 ±0.2 


0.052 ±0.005 


0.064 ±0.001 


Total R 


7.9 ±0.2 


1.41 ±0.05 


1.45 ±0.01 



a Total mass in combined system including the mass within both 
apertures, as well as outsi de the apertures in the image in which 
our S/N criteria (see 32.51) was met. 



in NGC 5195. Both galaxies reach stellar mass surface 
densities of over 2000 M pc -2 in the nuclear regions, 
then decrease radially to values of ~5OM pc -2 at the 
edge of the apertures (by definition) . The central surface 
density (in a 30" square region) is only slightly higher in 
the NGC 5195 at a value of 3490±900 M pc~ 2 compared 
to 2980 ± 450 M pc~ 2 in NGC 5194. 

The dust mass density in the central region of 
NGC 5194 peaks just above 2 M pc -2 in the central re- 
gion, decreasing to 0.75 M pc -2 along the spiral arms 
and then a constant 0.25 M pc -2 throughout the rest of 
the system. The dust mass density is comparably lower 
in NGC 5195 at values near 0.25 M pc -2 , consistent 
with being an extension of the spiral arm of NGC 5194. 

We provide the total stellar, gas and dust masses in 
each aperture in Table [2j We also measure the total mass 
in the combined system, summing up the mass in each 
pixel over the entire system. A small percentage (< 5%) 
of the galaxy is missing from these 'total' values as we 
only added u p pi xels which satisfied our signal-to-noise 
criteria (see £2.5) and low surface brightness regions of 
the galaxy are not accounted for. About ~ 80% of the 
gas and dust mass is found in the spiral NGC 5194, with 
the remaining found outside the aperture. Less than 1% 
of the gas and 3% of the dust of the total system is lo- 
cated within the aperture defining NGC 5195. The stellar 
mass is distributed more evenly between the two galax- 
ies with NGC 5194 being twice as massive as the spheroid 
NGC 5195. The bulk of the stellar mass is found within 
our defined apertures due to relatively steeper mass pro- 
files. Only about 10% of the stellar mass is located out- 
side the two apertures. 

4.1.2. Dust-to- stellar mass ratio 

Recall from Figures [5] and [9] which show the stellar 
and dust mass surface density maps of the system, that 
the dust mass is slightly more structured than the stel- 
lar mass map. In the stellar mass map, there are some 
enhancements in stellar mass along the spiral arms of 
NGC 5194, but they are not as pronounced as the spiral 
arms in the £?-band image (which can be found in the 
appendix) for example or those found in the dust mass 
map. Both the £?-band and dust mass surface density 
map are enhanced by the younger stellar populations. 
Because of this, the dust-to-stellar mass ratio, shown in 
Figure [l4j is a factor of two larger along the arms. 

We find an average dust-to-stellar mass ratio of 
M dust /M. = (3.1 ± 1.0) x 10- 3 (or log(M dust /M*) = 
-2.5±0.2) in NGC 5194. The dust-to-stellar mass ratio is 
an order of magnitude less in NGC 5195 at M^st/M* = 
(3.6 ± 2^ x 10- 4 (or log(M dust /M*) = -3.5 ± 0.3). In 
Figure [l5j the radial profiles of the dust-to-stellar mass 
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Fig. 13. — Radial mass profile of the stars (black dots), gas (blue 
triangles), and dust (red squares) in NGC 5194 (left) and NGC 5195 
(right). Each data point represents a 30" wide pixel box encom- 
passing one FWHM of the convolved beam size. Error bars are 
smaller than the symbol size if they are not visible. 



ratio is given for both galaxies. 

The dust-to-stellar mass ratio of NGC 5194 is consis- 
tent with results of other spiral galaxies in th e Herschel 



Reference Survey (HRS; |Boselli et al. 2Q1Q| ) w hich are 
found to ha ve log(M dust /AZ*) = -2.59 ± 0.03 flCortese 
et al.| |2Q12[ ) . Both galaxies are consistent with the 
mass ratio found for a large range of galaxy types, 



2.95 ±0.68 (|Skibba et al |2QlT]) 
vey (Kennicutt et al.||2U"ll|), g 



from 
a Her- 



log(M dust /M*) 
the KINGFISH survey 

schel extension of the SllNGS survey ( |Kennicutt et al.| 
|2QQ3[ ). Furthermore, the late- type NGC5 194 is consist 
tent w ith late-types in KINGFISH which Ski bba et al.| 
( 2QTTt found to be log(M dust /M*) = -2.70 ± 0.53 and 
NGC 5195 is in line with their findings for the early- type 
galaxies in the survey with \og(M^ ust / M^) = —3.77 ± 
0.83. Early- type galaxies in the volume based HRS show 
the spheroid population has a lower dust-t o-stellar mass 
ratio with log(M dust /M*) = -4.1 ± 0.1 (ISmith et al. 



2012[ ) compared to NGC 5195, although a couple ot ob- 



jects do have ratios as high as log(M dust /M*) = —3.5. 
For clarity, all the comparisons listed here have stellar 
and dust masses that w ere derived using a similar IMF 
(either a |Chabrier||2003| or |Kroupa||2001| IMF) and dust 
opacity and do not sutler from systematic offsets in the 
derived masses. 
In comparison to larger volume-based inf rared surveys, 



recen t results from the H- ATLAS survey (Bourne et al 
2012) reveal that for blue galaxies there is a trend tor 
the dust-to-stellar mass ratio to decrease with stellar 
mass. They adopt a higher opacity value than ours with 
^(A=250/im) = 0-89 m 2 kg _1 compared to our choice of 
^(A=250/im) = 0.398 m 2 kg -1 . For comparison, their dust 
masses and dust-to-stellar mass ratios need to be shifted 
up by 0.35 dex. Stellar m asses should be roughly consis- 
tent as their chosen IMF ( |Chabrier|2003[ ) leads to stellar 
masses only .03 dex lower than our chosen IMF from 



Kroupa| ( |200l| ). Accounting for this, NGC 5194's dust-to- 
stellar mass ratio is consistent with the blue galaxies from 
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Fig. 14. — A dust-to-stellar mass ratio map with the stellar mass 
density given in contours at the levels of 50, 100, 200, 500, 1000, 
1500, 2000, 2500, 3000 M pc~ 2 . NGC 5195's dust-to-stellar mass 
ratio is an order of magnitude smaller than NGC 5194 and increases 
with radius as the stellar mass decreases. 



the H-ATLAS survey (with log(M dust /M*) - -2.8), 
while NGC 5195 is slightly lower, consistent with the 
green galaxies in their sample, which it would belong 
to. 



4.1.3. Gas-to-dust mass ratio 



The map of the gas-to-dust ratio is shown in Figure [16 
along with the ID radial profile for NGC 5194. The ratio 
is fairly constant across NGC 5194 and appears to slightly 
decrease radially from a value of 105 ± 21 in the inner 
60" to a value of 71 ± 27 at a distance of 13 kpc from 
the galaxy's center. Recalling the metallicity map of the 
system in Figure [7] and metallicity gradient results from 
et al.| p( 



Moustakas et al. (2010), the metallicity drops by a factor 
1 1 



ot ^2 from t he central region to our aperture edge. Re- 



cent work by Ma grini et al.| ( |2011[ ) and |Foyle et al | ( |2012 ) 
show that a metallicity dependent Xqo conversion tac- 
tor is needed for regions dominated by molecular gas, 
which is the case for NGC 5194. Because of the radially 
decreasing metallicity, such a conversion would cause the 
gas-to-dust ratio to increase at larger radii leading it to 
be even more constant across the galaxy as was see n for 
M83, also observed in the VNGS (|Foyle et al.||2012[ ). 

In most cases, when a gradient is found in the gas-to- 
dust ratio, it often goes in the other direction, increas- 
ing with radial distance and is usually attributed to the 
galaxy containing a larger extended gas di sk relative to 
its dus t disk. This was seen in NGC 2403 ( Bendo et al. 
| 2010b[) and in the s ample of SI NGS galaxies ([Mu hoz- 
Mat eos et al.||2009j ). However, |Muhoz-Mateo? 



et al. 
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Fig. 15. — Radial profile of the dust-to-stellar mass ratios in 
NGC 5194 (purple) and NGC 5195 (cyan). The ratio is relatively 
constant in NGC 5194 but increases with radius in NGC 5195 pri- 
marily d ue to the decreasing stellar mass density as demonstrated 
in Figure [14] 



(|2009|) actually find a fairly constant gas-to-dust ratio 
tor NGC 5194 in agreement with our findings. 

4.2. Heating 
4.2.1. Decoupling of dust mass and heating 

There is a subtle decoupling between the distribution 
of the radiation intensities across the system and the dust 
mass surface density. Although both show spiral struc- 
ture, when we superimpose them together as shown in 
the rightmost plot of Figure [TtJ it appears that there is 
a slight offset between the peaks of each quantity. This 
has been s een already in the nearby spiral M83 (Foyle 
et al.|[2012 ), also observed as part of the VNGS, as well 
asm NGC 4501 an d NGC 4567/8 in th e Virg o Cluster 



([Smith et aLl|201Q[). As in | Foyle et al.| ([2012), we find 



the peaks in the 1SRF (or dust temperature values) lie 
downstream of the spiral arm distribution of the dust, 
assuming a trailing arm rotation of NGC 5194 (counter- 
clockwise). 

There are two more obvious discrepancies between the 
U m in map and the dust mass map in the central regions 
of both galaxies. In the center of NGC 5194, there is 
an obvious increase in the dust mass density, relative to 
both the large scale mass density as well as the mass den- 
sity in the spiral arms. On the other hand, the ISRF is 
enhanced but only at the same level as the spiral arms. 
Even though both the stellar, dust and gas mass densities 
are all increasing, the local ISRF value stays moderately 
low. This is unlike the central region of NGC 5195 which 
has extremely high ISRF values (with peak dust temper- 
atures of 30 K), but no relative increase in the dust mass. 
Although there is less dust in the vicinity of NGC 5195, 
the dust that is there is quite warm. This is in a greement 
wit h mea surements of the gas temperature by |Kohno et 
aD ([2002) who found that low HCN-CO ratios suggest 
that the molecular gas in NGC 5195 is quite warm and 
unable to collapse and form stars. 



4.2.2. Ionized gas and the ISRF 

In the left panel of Figure [lT| we show the narrow- 
band Ha image with stellar emission subtracted. This 
was done by calculating the stellar continuum transmit- 
ted through the Ha narrowband filter curve using the 
best-fit SED model without nebular emission included. 
The strength of the ISRF, as parametrized by U : 



is 



10]). Qualitatively, 
Ty correlated in the 



shown in contours (and also in Fig. 
the Ha emission and ISRF are spatial 
disk of NGC 5194. They are both enhanced along the 
spiral arms and in the central bulge region, although Ha 
is more peaked than the ISRF in the nucleus perhaps due 
to AGN emission. In contrast, in NGC 5195, there is no 
detected star formation and hence no Ha emission. We 



note howeve r , that in optical spectroscopy from |Mous- 
takas et al. (2010), a modest amount of Ha emission 
is detected in the nuclear region at a ~2a level that ap- 
pears to be blue shifted relative to Ha absorption in stel- 
lar photospheres. Surprisingly, this is where the ISRF is 
the brightest in the entire map. In Figure [18| we show a 
ID comparison between the ionized gas intensity and the 
strength of the ISRF determined from the dust emission 
model fitting. While the ionized gas intensity correlates 
with the ISRF in NGC 5194, it clearly follows another 
trend in NGC 5195 where high ISRF values correspond 
to no Ha emission. 

4.2.3. Dust heating in NGC 5195 

One of the most interesting findings of this analysis is 
that the ISRF of the post-starburst spheroid NGC 5195 
is up to 20-25 times the characteristic amplitude of the 
ISRF in the Milky Way despite its lack of star formation. 
Often high dust temperatures are associated with star- 
bursting and actively star forming galaxies. This is not 
the case for NGC 5195. The galaxy has a high (29 ± 3 K) 
characteristic dust temperature in the nuclear region and 
an average temperature of 23 ± 3 K in the galaxy's aper- 
ture. To provide some context, recen t results from the 
H- ATLAS survey ( |Bourne et al.||2012[ ) indicate that, for 
blue galaxies in their sample, dust temperatures increase 
with stellar mass up to a stellar mass of ~ 6 x 10 10 M 
reaching a peak dust temperature of ~ 25 — 30 K and 
then anti-correlate at higher stellar masses (see their Fig- 
ure 14). NGC 5194, a blue galaxy according to their 
criteria, is consistent with this trend as is NGC 5195. 
However, NGC 5195 is not a blue galaxy, but is about 
0.4 mag redder than NGC 5194 and would thus fall into 
their green galaxy color bin. In this colour bin, a few 
galaxies exhibit high temperatures similar to NGC 5195, 
but their temperatures are poorly determined likely be- 
cause the dust model used to estimate the temperature 
was a poor fit to the observations. The typical scenario 
suggests the redder the galaxy, the cooler the dust tem- 
perature, but there are hints, particularly from the large 
scatter and uncertainty in constraining the temperature, 
that the transition from blue-to-red galaxy is more com- 
plex in the context of dust heating. 

As a post-starburst galaxy, NGC 5195 is in a transi- 
tion from starburst to q uiescent. The lack of ionized gas 
(Thronson et al. 1991) indicates no recent star forma- 
tion is ongoing but did peak during its close encounter 
with NGC 5194 ~ 400 Myr ago. Other signposts that it 
is in a transition are evident by its low dust-to-stellar 
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and gas-to-stellar mass ratios. Its gas fraction at ~2% 
is very low and we show, in Figure [l9j the gas that is 
in the galaxy appears to mostly be below the empirical 
threshold o bserved in the THINGS galaxy survey ( |Bigiel| 
et al.|2QQ8[ ) at which atomic gas transitions to molecular 
gas, a requirement for star formation. In this figure we 
show the stellar mass density as a function of gas mass 
density in the top panels, with NGC 5194 in the top- left 
and NGC 5195 in the top-right. The axes are kept on the 
same scale for each panel to show the different density 
parameter space in each galaxy. For a given stellar den- 
sity, the gas density is considerably lower in NGC 5195 
and there is a large scatter between stellar mass density 
and the total gas mass density (purple points). On the 
other hand, NGC 5194 shows a tighter relation between 
the total gas (purple) and stellar mass densities and there 
appears to be a threshold gas density of ~lOM pc -2 
(marked by vertical dashed lines) at which the atomic gas 
(open red squares) is converted from neutral to molec- 
ular (blue plus signs). This threshold was noticed in 
the Kennicutt- Sch midt (K-S) relations drawn from the 
THINGS survey in |Bigiel et H] fl2008p . In their plots 
they used star formation density on the y-axis and found 
the same transition from atomic to molecular gas at a 
fixed gas density of 9M pc -2 . In the bottom panel, we 
provide plots analogous to these K-S plots by plotting 
the intensity of Ha emission (a proxy for star formation 
rate density) in our narrowband images after subtracting 
stellar continuum based on our SED fitting. We do not 
correct for dust attenuation or conve rt the intensity to 



i 2 



2008) 



SFR density. The reader is referred to pigiel et al 
who already performed this analysis for NGC 5l0 
NGC 5195, no star formation is found in the galaxy so 
no Kennicutt-Schmidt relation can be observed in this 
galaxy. 

So why is the dust temperature so high? We attempt 
to shed some light on the mechanism leading to the high 
temperature observed in this galaxy. Because this galaxy 
represents a nearby transitional galaxy, studying its re- 
solved multi-wavelength properties can help isolate the 
source of heating. In the following few sections we will 



explore some of the possible heating mechanisms: (1) 
heating from star formation, (2) heating from the evolved 
stellar population and (3) heating by its active galactic 
nucleus. 

Heating by ongoing star formation: It is clear that the 
ionizing radiation of young massive stars is not heating 
the dust as evident by its lack of ionized gas shown in 
Figure 18 The pixels in NGC 5195 are all upper limits in 
Ha so no ionized gas is detected except for a few regions 
near the edge of the galaxy, at which the emission from 
NGC 5195 begins to blend with the spiral and the local 
stellar density is lower. 

Heating by t he ev olve d population: Rece nt studies by 
Bendo et al Jfl2012| ) and |Boselli et al.| ( |2012| ) indicate that 
evolved stellar population heating plays a strong role in 
dust heating, although this is most important for lower 
temperature dust heating at longer infrared wavelength 
beyond 160 fim. In Figure [20j we plot the amplitude of 
the ISRF as a function of stellar, total gas and dust mass 
density. The ISRF does not correlate with either the 
stellar, dust or gas mass density. In addition, the bulge 
in the spiral reaches similarly high stellar mass densi- 
ties (~ 3 x 10 4 M o pc -2 ) but has a much lower ISRF 
value compared to NGC 5195, so the enhanced stellar 
mass density cannot be the primary explanation for the 
enhanced ISRF observed in NGC 5195. 

Heating by its active galactic nucle us: NGC 5195 has a 
comp act radio source in the nucleus ( |van der Hulst et aL 
| 1988[) and has spectral indications of harboring an AGN 
( jlVloustakas et al.|2010 ). In a spectrum presented in cov- 
ering the nuclear region, a small amount of Ha emission 
is seen blueshifted relative to Ha seen in absorption, in- 
dicative of an outflow of ionized gas (Moustakas et al. 
2010). In our observations, the ISRF is strongly peaked 
in the nuclear region and is suggestive that the AGN 
plays a role in heating the dust. 

However, it is not completely clear that this is the dom- 
inant mechanism heating the dust and/or shutting down 
star formation. Mid-infrared observ ations with ISOCAM 
show a lack of AGN emission lines ( |Boulade et al.|[l 996) 
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Fig. 17. — Left: Map of the stellar subtracted Ha intensity with U m i n shown in contours. The highest ISRF is found in the post-starburst 
galaxy, NGC 5195, where a clear lack of Ha emission is found in the nuclear region. The ISRF is enhanced along the spiral arms, as is Ha 
emission but does not show an increase in the nuclear regions, in contrast to the Ha emission which is strongly peaked in the bulge. Right: 
Mass of the dust mass surface density in Mq pc -2 with contours of the amplitude of the ISRF of the ISM, quantified by U m i n from DL07 
models. While both are enhanced along the spiral arms, there are some noticeable anti-correlations, specifically in the nuclear regions of 
both galaxies. 
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Fig. 18. — The Ha emission is moderately correlated with the 
ISFR in NGC 5194 (purple points), while NGC 5195 (cyan triangles 
and upper limits) has minimal to zero Ha emission despite having 
high ISRF values. 



and rule out LINER activity as the dominant source of 
mid-infrared emission. Furthermore, NGC 5194 also har- 
bors an AGN and bas ed on equivalent width s and line in- 
tensities presented in Moustakas et al. (2010), the ionized 
emission associated with the AGN is 5-10 times brighter 
in NGC 5194. We caution that this comparison is compli- 
cated by unknown dust a ttenuation due t o strong Balmer 
absorption in NGC 5195 ( [Ho et al.|1995| ). Unlike the nu- 



clear region of NGC 5195, the ISRF is only slightly en- 
hanced in the nuclear region in NGC 5194 (see left panel 
in Figure 17), and at the same intensity as is seen in the 
spiral arms. 

Unfortunately none of these mechanisms stand out as 
the clear heating source. Although evolved population 
heating appears to be the most likely, with the other two 
almost ruled out. The stellar mass density is similar in 
both galaxies so have to conclude that the relatively lower 
gas and dust fractions in NGC 5195 result in less efficient 
gas and dust cooling, leading to a relatively higher ISRF. 

4.2.4. PAH destruction 

Looking back on Figure [ToJ the PAH fraction decreases 
slightly to ~ 4% in the nuclear region of NGC 5194 due 
to a low-lu minosity active galactic nucleus (Moustakas 
et al. 2Q1Q|). This is not too surprising. M any studies 
report ( |Uenzel et aL]|1998 Dale et ar]|2QQ6[ ) the destruc- 
tion PAH molecules aro und AGN a nd withi n hard radi- 
ation fields with PDRs jBos elli et aLl|2QQ4| |Calzetti et 

On the other hand, the de- 



aLl[2QQ5l |Bendo et al.| 2008 ) . 
crease m PAH traction in N 



GC 5195 is interesting. The 
PAH fraction traces the stellar mass surface density and 
decreases as a gradient across NGC 5195. One can see 
in the bottom-right SED of NGC 5195 in Figure [4] that 
PAH emission is easily detected and is many times the 
stellar emission at 5.8 and 8.0 jam. So the decrease in the 
mass fraction of PAHs relative to the total dust mass is 
not a detection issue. N GC 5195 also harbours an AGN 
( |Moustakas et al.|2010 ) leading to a lower PAH fraction. 
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Fig. 19. — Stellar mass density versus atomic (red squares), molecular (blue plus signs) and total (purple dots) gas mass surface density 
for NGC 5194 in the top-left and for NGC 5195 in the top-right. Stellar subtracted Ha intensity as a function gas mass surface density is 
shown in the bottom panels. In each panel, we note the surface density threshold of 10 Mq pc - 2 at which atomic gas converts to molecular 
gas. Note that some regions in NGC 5195 contain molecular gas at densities lower than this threshold. 

But unlike NGC 5194, where the depression is concen- 
trated in the nuclear region, the decrease in PAH frac- 
tion across NGC 5195 is a gradient. This may suggest 
NGC 5195 is accreting gas, dust and PAH material from 
the nearby spiral arm, but as the PAHs approach the 
AGN, the fraction decreases. 



formation is completely suppressed ( |Kennicutt|1989| ) be- 
cause the molecular gas does not become dense enough 
to collapse into cold, dense molecular clouds required for 
star formation. 
This scenario is supporte d by observations of the 



4.3. An example of morphological quenching 

Unlike NGC 5194, the molecular gas in NGC 5195 is 
found at densities lower than ~lOM pc -2 . If this 
threshold marks the density at which atomic gas be- 
comes molecular, then it is likely that the molecular gas 
has come from NGC 5194 via accretion rather than be- 
ing formed in situ. In addition to this, the molecular 
gas that does exist is unable to form stars because of the 
high stellar density (and thus its gravitational potential). 
The critical gas density required for the gas to cool and 
collapse is not satisfied. Mo delled as a simple single- 
fluid T oomre stability model ( |Toomre||1964 ), Kennicutt 



molecular gas in NGC 5195 b yjKohno et~al| ( |2002| ). They 
measured both CO (1-0) and H (JIN (1-0) in the galaxy and 
found an HCN-to-CO intensity rati o 5-15 times smaller 



than that seen in starburst regions ( |Helfer fc Blitz||1993 i 
Kohno et al.|2 002) an d 5 times smaller than that usually 



found in spiral disks (Heifer & Blitz||1995 Kohno et al. 



(1989) show that at higher stellar densities, the molec 
ular gas must reach a higher critical density in order to 
cool and collapse. This is because localized shear motions 
from the gravitational potential are greater than the self- 
gravitational potential of the molecular cloud thus mak- 
ing the gas stable. In regions where the average surface 
density of gas is less than 5 times the critical density, star 



1999), indicative of a lack ol dense molecular cores 

Half a billion years ago, when the two galaxies were 
closer together, tidal forces would have disrupted such 
high shear velocities and likely involved higher rates of 
gas and dust accretion so both the density required for 
clumps to form would have been lower and the gas sup- 
ply available would have been more dense, allowing star 
formation to occur. But as the galaxies drew apart, the 
tidal forces were reduced and the dense molecular gas 
was exhausted as it was converted into stars during the 
galaxy's starburst episode. 

This process, recently ter med "morphological quench- 
ing" by Martig et al. (2009), occurs naturally within hy- 
dro dynaim^aTimiuTiitions of galaxy evolution and c an be 
extremely efficient under the right circumstances. Mar- 
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Fig. 20. — The amplitude of the ISRF as a function of stellar mass (left), total gas mass (middle), and dust (right) for NGC 5194 (purple 
dots) and IN GC 5195 (cyan triangles). There is a modest trend lor higher values to be lound at higher stellar densities. 



tig et aL] ( |2QQ9| ) showed that gaseous molecular disks, 
even those up to several 10 9 M , are unable to form stars 
if they are embedded within a sufficiently large grav- 
itational potential such as a bulge dominated galaxy. 
Star formation will cease even if gas accretion contin- 
ues. The authors claim that this mode of star formation 
quenching is important for galaxies in dark matter ha- 
los less massive than ~ 10 12 M and thus an important 
mechanism for converting L* and lower mass galaxies 
from the 'blue cloud' to the 'red-sequence' and regulat- 
ing galaxy growth, particularly because other quenching 
mechanisms such as AGN feedback or gas stripping are 
only effective in more massive halos. NGC 5195 provides 
a great nearby laboratory to further understand and val- 
idate this star formation quenching mechanism. 

5. CONCLUSIONS 

We have fit optical, near-, mid- and far- infrared im- 
ages of the Whirlpool galaxy system to a large library of 
stellar and dust spectral energy distribution models. By 
performing 100 Monte Carlo simulations of the photom- 
etry fitted to the models in each pixel, we have obtained 
maps of output parameters from st ellar models generated 
by th e SPS code PEGASE.2 flFioc fc Rocca-Volmerange 
|1997| ) and dust m odel parameters gen erated by the dust 
emission model of |Draine fc Li| (2007). Our images were 
matched to a common plate scale and spatial resolution 
defined by a gaussian PSF with a FWHM=28 // resolv- 
ing the galaxy on physical scales of ~ 1 kpc. Here we 
summarize the main conclusions of this work: 

1. The SPS model fits reveal a burst of star formation 
occurred in both galaxies roughly 340-500 Myr ago, 
consistent with dynamica l models of the interaction 



history of the gala xies (Salo & Laurikainen 2000 
Dobbs et al.|201Q ) as well as colour-magnitude di - 



agrams of individual stars ( Tikhonov et al.||2009| ). 
There is little spatial dependence of the burst age 
and other SFH output parameters, suggesting a 
resolved SED approach is not necessary at this 
spatial resolution, corresponding to physical scales 
of ~lkpc. Metallicity and dust attenuation show 
some radial dependence but were not very well con- 
strained by our method. 

2. The dust-to-stellar mass ratio in NGC 5194 is 



log(M dust /M*) = -2.5 ± 0.2, while NGC 5195 has 
a magnitude lower mass ratio at NGC 5195 at 
log(Md us t/M*) = — 3.5±0.3. Both galaxies are con- 
sistent, according to galaxy t ype, with other galax- 



ies observed with Herschel (iCo rtes e et al. 2012 
Smith et al.|2012| |Bourne et aTpOT?| ~ 



3. High radiation field values concentrated in the nu- 
clear region of NGC 5195 provide evidence that the 
dust is located within the galaxy, not in its fore- 
ground. Despite a high ISRF (often an indication 
of star formation), there is no current star forma- 
tion in the central region of NGC 5195 as indicated 
by a lack of Ha emission and l ack of dense molec- 
ular clouds flKohno et aT1|2002 ). 



4. Morphological quenching ( Martig et al.||2009| ) is a 
plausible scenario for the suppression of star lorma- 
tion in NGC 5195. With its higher gravitational po- 
tential, the gas density threshold required for star 
formation is higher than the available low-density 
gas and is thus unable to form stars. 
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APPENDIX 

For reference, we show the image processing of the archival data used for our analysis. As described in £p} a ll images 
were matched to a common spatial resolution (FWHM=28 // ) and plate scale (10" pix -1 ). Fig ures [2l| fc |22| show the 
optical images taken from the ancillary catalogue of the SINGS survey ( Kennicutt et al.||2003|) . TheTeft image is the 
original image of the galaxy at native resolution prior to any post-processing. The middle panel shows the galaxy after 
foreground stellar emission is replaced by emission in surrounding pixels. Following this, we are able to match the 
spatial resolution to the far-infrared images using convolution kernels as described in ^2] The far right image shows a 
normalized residual map which we quantify as in Equation [2] 

Most of the optical images are well fit by the modes with the residual images revealing that the models fit the data 
to within 10%, and thus within their uncertainties. The only image that is not well fit is the Ha image. This image 
was not used in the fitting routine but since the SEDs include nebular emission, ideally, it should be predicted by 
the SED models. The predicted nebular emission is reasonable in NGC5195, which has no star formation, but in 
NGC 5194, the amount of Ha emission is underestimated. This is primarily because our chosen SFH is too simplistic 
and only modelled the very old population and the burst population from ~ 400 Myr ago. A more complex SFH, such 
as a declining or constant star forming burst, is needed to properly describe the recent star formation occurring in 
NGC 5194. 

The near-infrared images from the 2MASS Large Galaxy Atlass (|Jarrett et al. [2003) are shown in Figure 23 The 



models lead to small residual images, indicating a good fit to the NIK images. On the other hand, mid-infrarecTTrnages 
are not as well fit as shown in Figure [24] This is primarily because the mid-infrared images at 3.6 and 4.5 /am are not 
used to fit the SEDs and the emission is being extrapolated from models that have fit on one side of the spectrum the 
stars and on the other side the dust and PAH emission. Both models combined tend to overestimate the emission at 
3.6 and 4.5 /am. More analysis is needed to determine whether the stellar or PAH emission is being overestimated or 
both. The 5.6 and 8.0 /am images are used to fit the DL07 dust emission SED model. In the nuclear regions of both 
galaxies, the models are a good representation of the emission, but outside of these regions, the models underestimate 
the emission. Recalling from the PAH mass fraction map on the right in Figure [10] that in the regions outside of 
the nuclear region or edges of the galaxies, the PAH mass fraction reaches the highest value permitted by the DL07 
models. This is likely the reason why the 5.6 and 8.0 /am are underestimated by the models. A higher PAH fraction 
is needed to properly describe the observations. At longer wavelengths, just like our Herschel observations shown in 
Figure |2j the 24 /am image is accurately described by the DL07 models. 
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Fig. 21. Optical images of NGC 5194/5195 from the SINGS survey ( |Kennicutt et al.|2003| > from top to bottom: B,V and Ha. On the 
left we show the original image at its native resolution in each band and then the processed images after convolution with the beam size of 
each filter in the bottom right of each panel. Both left and middle panels are in intensity units of MJy sr — 1 , whose range is indicated by the 
gradient on the right of the middle panel. The right column shows the residual image, defined in Eq. [2] between the observed image after 
it is PSF matched (middle) and a synthetic galaxy image (not shown) generated by the best fitting SED in our analysis. In the bottom 
row, the observed and processed Ha narrow band image is shown (and includes stellar continuum emission in this band as there was no 
need to stellar continuum subtract in our analysis). The image was not include in our SED fitting process, but we do show that overall the 
models - which include nebular emission- underestimate the amount of nebular emission by more than 20%, except in NGC 5195, where 
the models accurately predict the amount of emission (mostly stellar continuum emission in this case) in the Ha bandpass. 
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Fig. 22. Optical images of NGC 5194/5195 from the SINGS survey ( |Kennicutt et al.||2003| from top to bottom: R and I. As in 
Figure pi] the left panel is the original image from the archive and the processed image in the middle, both in MJy sr _1 . The right shows 
the residual image between the models and observations.. 
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Fig. 23. Near-infrared JHK images of NGC 5194/5195 from the 2MASS Large Galaxy Atlas jJarrett et al.|2003| >. As in Figure [21) the 
left panel is the original image from the archive and the processed image in the middle, both in MJy sr _i . The right shows the residual 
image between the models and observations. 
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Fig. 24. Mid-infrared images of NGC 5194/5195 from Spitzer as part of the SINGS survey flKennicutt et al.||2003| ). As in Figure |2l] 
the left panel is the original image from the archive and the processed image in the middle, both in MJy sr _i . The right shows the residual 
image between the models and observations. The top two images were not used in either fit. 



